********************************************************************************
* 						 CHEMICAL WEAPONS ANALYSIS							   *
********************************************************************************

clear all
set more off
set scheme plotplainblind
macro drop _all
est drop _all
set matsize 800
set seed 8675309

** Set Working Directory

if c(username) == "christopherblair"{
global dir "~/Desktop/JOP Replication"
global raw "${dir}/Raw Files"
global code "${dir}/Code"
global result "${dir}/Results"									
}

else if c(username) == "youruser"{
global dir "~/Desktop/JOP Replication"
global raw "${dir}/Raw Files"
global code "${dir}/Code"
global result "${dir}/Results"	
}

cd "$raw"

********************************************************************************

use "${dir}/ams175.dta", clear

sort ballotnum

global demographic "age school monthsoversea rankgrade"
global covariates "ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact"

********************************************************************************
*								GRAPHICAL RESULTS							   *
********************************************************************************

collapse (mean) meanuse= use (sd) sduse=use (count) n=use, by(numchamber)
generate hiuse = meanuse + invttail(n-1,0.025)*(sduse / sqrt(n))
generate lowuse = meanuse - invttail(n-1,0.025)*(sduse / sqrt(n))

twoway (bar meanuse numchamber, fcolor(white) lcolor(black)) (rcap hiuse lowuse numchamber, lcolor(black)), ylabel(.7(.05)1) ymtick(.7(.01)1) yline(1, lcolor(gs10) lpatt(dot)) yline(.9069401, lcolor(cranberry) lpatt(shortdash)) ytitle("Average Support for Using" "Chemical Weapons Against Japan") xtitle("Times Training w/ Gas Exposure") legend(off) text(.925 -.22 "0.933", size(small)) text(.897 .82 "0.907", size(small)) text(.8875 1.82 "0.897", size(small)) text(.865 2.82 "0.871", size(small))
graph export "${result}/mean_chamber_ci.png", replace

use "${dir}/ams175.dta", clear

sort ballotnum

global demographic "age school monthsoversea rankgrade"
global covariates "ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact"

collapse (mean) meanuse= use (sd) sduse=use (count) n=use, by(nummask)
generate hiuse = meanuse + invttail(n-1,0.025)*(sduse / sqrt(n))
generate lowuse = meanuse - invttail(n-1,0.025)*(sduse / sqrt(n))

twoway (bar meanuse nummask, fcolor(white) lcolor(black)) (rcap hiuse lowuse nummask, lcolor(black)), ylabel(.7(.05)1) ymtick(.7(.01)1) yline(1, lcolor(gs10) lpatt(dot)) yline(.9069401, lcolor(cranberry) lpatt(shortdash)) ytitle("Average Support for Using" "Chemical Weapons Against Japan") xtitle("Times Training in Gas Masks") legend(off) text(.8875 -.22 "0.899", size(small)) text(.818 .82 "0.828", size(small)) text(.968 1.82 "0.979", size(small)) text(.897 2.82 "0.906", size(small)) text(.91 3.81 "0.919", size(small))
graph export "${result}/mean_mask_ci.png", replace


********************************************************************************
*								SUMMARY STATISTICS							   *
********************************************************************************

use "${dir}/ams175.dta", clear

sort ballotnum

global demographic "age school monthsoversea rankgrade"
global covariates "ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact"

foreach x of var school age monthsoversea rankgrade {
tab `x', gen(`x')
}

global school "school2 school3 school4 school5 school6"
global age "age2 age3 age4 age5 age6 age7 age8 age9 age10 age11 age12"
global rank "rankgrade2 rankgrade3 rankgrade4 rankgrade5 rankgrade6"
global months "monthsoversea2 monthsoversea3 monthsoversea4 monthsoversea5 monthsoversea6 monthsoversea7 monthsoversea8 monthsoversea9 monthsoversea10"
global dv "use firstuse seconduse"

estpost sum $dv numchamber nummask $school $age $rank $months $covariates
esttab using "${result}/CW_SumStat.tex", style(tex) cells("count(pattern(1 1 0) fmt(0)) mean(pattern(1 1 0) fmt(3)) sd(pattern(1 1 0) fmt(3)) min(pattern(1 1 0) fmt(3)) max(pattern(1 1 0) fmt(3))") varlabels(use "& \\ \textsc{Dependent Variables} \\ \hspace{3mm}Support Using Chemical Weapons Against Japan" firstuse "\hspace{3mm}Support First-Use Against Japan" seconduse "\hspace{3mm}Support Second-Use Against Japan" numchamber "& \\ \textsc{Independent Variables} \\ \hspace{3mm}\# of Times Exposed to Gas" nummask "\hspace{3mm}\# of Times Training w/ Gas Mask" school2 "& \\ \textsc{Control Variables} \\ \hspace{3mm}Schooling: <\nth{8} Grade" school3 "\hspace{3mm}Schooling: \nth{8} Grade" school4 "\hspace{3mm}Schooling: Some High School"  school5 "\hspace{3mm}Schooling: High School" school6 "\hspace{3mm}Schooling: College" age2 "\hspace{3mm}Age: $\leq$18" age3 "\hspace{3mm}Age: 19" age4 "\hspace{3mm}Age: 20" age5 "\hspace{3mm}Age: 21" age6 "\hspace{3mm}Age: 22" age7 "\hspace{3mm}Age: 23" age8 "\hspace{3mm}Age: 24" age9 "\hspace{3mm}Age: 25" age10 "\hspace{3mm}Age: 26-29" age11 "\hspace{3mm}Age: 30-34" age12 "\hspace{3mm}Age: 35+" monthsoversea2 "\hspace{3mm}Time Deployed: $\leq$3 Months" monthsoversea3 "\hspace{3mm}Time Deployed: 3-6 Months" monthsoversea4 "\hspace{3mm}Time Deployed: 6-9 Months" monthsoversea5 "\hspace{3mm}Time Deployed: 9-12 Months" monthsoversea6 "\hspace{3mm}Time Deployed: 12-18 Months" monthsoversea7 "\hspace{3mm}Time Deployed: 18-24 Months" monthsoversea8 "\hspace{3mm}Time Deployed: 24-30 Months" monthsoversea9 "\hspace{3mm}Time Deployed: 30-36 Months" monthsoversea10 "\hspace{3mm}Time Deployed: 36+ Months" rankgrade2 "\hspace{3mm}Rank: PRV./PFC." rankgrade3 "\hspace{3mm}Rank: CPL./TCH5." rankgrade4 "\hspace{3mm}Rank: SGT./TCH4." rankgrade5 "\hspace{3mm}Rank: SSGT./TCH3." rankgrade6 "\hspace{3mm}Rank: TSGT./MSGT./1SGT." ally_relations "\hspace{3mm}Postwar Foreign Policy Index" orientation_officers "\hspace{3mm}Officers' Leadership Index" infocenter "\hspace{3mm}Unit Has War Information Center" orient_meet "\hspace{3mm}Unit Holds Regular Orientation Meetings" war_interest "\hspace{3mm}Interest in War News Index" honolulu_contact "\hspace{3mm}Local Contact Index") nolines prehead(\begin{tabular}{l*{7}{c}} \hline & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) collabels("Observations" "Mean" "Std. Dev." "Min" "Max") nomti nonum noobs replace


********************************************************************************
*								DIFFERENCE-IN-MEANS							   *
********************************************************************************

eststo sum0: qui estpost summ $school $age $rank $months $covariates if chamber == 0
eststo sum1: qui estpost summ $school $age $rank $months $covariates if chamber == 1
eststo diff: qui estpost ttest $school $age $rank $months $covariates, by(chamber)

esttab sum0 sum1 diff using "${result}/CW_DiffMeans1.tex", style(tex) cells("mean(pattern(1 1 0) fmt(3)) b(star pattern(0 0 1) fmt(3))" "sd(pattern(1 1 0) par)") varlabels(school2 "& \\ \hspace{3mm}Schooling: <\nth{8} Grade" school3 "\hspace{3mm}Schooling: \nth{8} Grade" school4 "\hspace{3mm}Schooling: Some High School" school5 "\hspace{3mm}Schooling: High School" school6 "\hspace{3mm}Schooling: College" age2 "\hspace{3mm}Age: $\leq$18" age3 "\hspace{3mm}Age: 19" age4 "\hspace{3mm}Age: 20" age5 "\hspace{3mm}Age: 21" age6 "\hspace{3mm}Age: 22" age7 "\hspace{3mm}Age: 23" age8 "\hspace{3mm}Age: 24" age9 "\hspace{3mm}Age: 25" age10 "\hspace{3mm}Age: 26-29" age11 "\hspace{3mm}Age: 30-34" age12 "\hspace{3mm}Age: 35+" monthsoversea2 "\hspace{3mm}Time Deployed: $\leq$3 Months" monthsoversea3 "\hspace{3mm}Time Deployed: 3-6 Months" monthsoversea4 "\hspace{3mm}Time Deployed: 6-9 Months" monthsoversea5 "\hspace{3mm}Time Deployed: 9-12 Months" monthsoversea6 "\hspace{3mm}Time Deployed: 12-18 Months" monthsoversea7 "\hspace{3mm}Time Deployed: 18-24 Months" monthsoversea8 "\hspace{3mm}Time Deployed: 24-30 Months" monthsoversea9 "\hspace{3mm}Time Deployed: 30-36 Months" monthsoversea10 "\hspace{3mm}Time Deployed: 36+ Months" rankgrade2 "\hspace{3mm}Rank: PRV./PFC." rankgrade3 "\hspace{3mm}Rank: CPL./TCH5." rankgrade4 "\hspace{3mm}Rank: SGT./TCH4." rankgrade5 "\hspace{3mm}Rank: SSGT./TCH3." rankgrade6 "\hspace{3mm}Rank: TSGT./MSGT./1SGT." ally_relations "\hspace{3mm}Postwar Foreign Policy Index" orientation_officers "\hspace{3mm}Officers' Leadership Index" infocenter "\hspace{3mm}Unit Has War Information Center" orient_meet "\hspace{3mm}Unit Holds Regular Orientation Meetings" war_interest "\hspace{3mm}Interest in War News Index" honolulu_contact "\hspace{3mm}Local Contact Index") nolines prehead(\begin{tabular}{l*{5}{c}} \hline & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1 2) collabels(" " " " " ") mti("Not Exposed to Gas" "Exposed to Gas" "Difference-in-Means (Not-Exposed - Exposed)") nonum noobs replace

eststo clear

eststo sum0: qui estpost summ $school $age $rank $months $covariates if multiple == 0
eststo sum1: qui estpost summ $school $age $rank $months $covariates if multiple == 1
eststo diff: qui estpost ttest $school $age $rank $months $covariates, by(multiple)

esttab sum0 sum1 diff using "${result}/CW_DiffMeans2.tex", style(tex) cells("mean(pattern(1 1 0) fmt(3)) b(star pattern(0 0 1) fmt(3))" "sd(pattern(1 1 0) par)") varlabels(school2 "& \\ \hspace{3mm}Schooling: <\nth{8} Grade" school3 "\hspace{3mm}Schooling: \nth{8} Grade" school4 "\hspace{3mm}Schooling: Some High School" school5 "\hspace{3mm}Schooling: High School" school6 "\hspace{3mm}Schooling: College" age2 "\hspace{3mm}Age: $\leq$18" age3 "\hspace{3mm}Age: 19" age4 "\hspace{3mm}Age: 20" age5 "\hspace{3mm}Age: 21" age6 "\hspace{3mm}Age: 22" age7 "\hspace{3mm}Age: 23" age8 "\hspace{3mm}Age: 24" age9 "\hspace{3mm}Age: 25" age10 "\hspace{3mm}Age: 26-29" age11 "\hspace{3mm}Age: 30-34" age12 "\hspace{3mm}Age: 35+" monthsoversea2 "\hspace{3mm}Time Deployed: $\leq$3 Months" monthsoversea3 "\hspace{3mm}Time Deployed: 3-6 Months" monthsoversea4 "\hspace{3mm}Time Deployed: 6-9 Months" monthsoversea5 "\hspace{3mm}Time Deployed: 9-12 Months" monthsoversea6 "\hspace{3mm}Time Deployed: 12-18 Months" monthsoversea7 "\hspace{3mm}Time Deployed: 18-24 Months" monthsoversea8 "\hspace{3mm}Time Deployed: 24-30 Months" monthsoversea9 "\hspace{3mm}Time Deployed: 30-36 Months" monthsoversea10 "\hspace{3mm}Time Deployed: 36+ Months" rankgrade2 "\hspace{3mm}Rank: PRV./PFC." rankgrade3 "\hspace{3mm}Rank: CPL./TCH5." rankgrade4 "\hspace{3mm}Rank: SGT./TCH4." rankgrade5 "\hspace{3mm}Rank: SSGT./TCH3." rankgrade6 "\hspace{3mm}Rank: TSGT./MSGT./1SGT." ally_relations "\hspace{3mm}Postwar Foreign Policy Index" orientation_officers "\hspace{3mm}Officers' Leadership Index" infocenter "\hspace{3mm}Unit Has War Information Center" orient_meet "\hspace{3mm}Unit Holds Regular Orientation Meetings" war_interest "\hspace{3mm}Interest in War News Index" honolulu_contact "\hspace{3mm}Local Contact Index") nolines prehead(\begin{tabular}{l*{5}{c}} \hline & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1 2) collabels(" " " " " ") mti("Not Exposed to Gas" "Exposed to Gas" "Difference-in-Means (Not-Exposed - Exposed)") nonum noobs replace

drop school1-rankgrade6

********************************************************************************
*									MAIN ANALYSES							   *
********************************************************************************

eststo clear
	
eststo: reghdfe use numchamber nummask, cluster(base) noabs
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber nummask, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber nummask ally_relations, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, vce(robust) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local Cluster "No"
	
esttab using "${result}/CW_Main.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic cov Cluster Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}Airbase Clustered SEs" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{17}{c}} \hline & \\ & & \multicolumn{15}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-17} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1 2 3 4 5 6 7 8) nomtit replace

eststo clear

********************************************************************************
*									OSTERS BOUNDS							   *
********************************************************************************

eststo clear
	
qui reg use numchamber nummask, cluster(base)
psacalc delta numchamber, rmax(.00676)

qui reg use numchamber nummask i.($demographic), cluster(base)
psacalc delta numchamber, rmax(.0637)

qui reg use numchamber nummask i.($demographic) ally_relations, cluster(base)
psacalc delta numchamber, rmax(.08788)

qui reg use numchamber nummask i.($demographic) ally_relations orientation_officers, cluster(base)
psacalc delta numchamber, rmax(.0897)

qui reg use numchamber nummask i.($demographic) ally_relations orientation_officers infocenter orient_meet war_interest, cluster(base)
psacalc delta numchamber, rmax(.10439)

qui reg use numchamber nummask i.($demographic) ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base)
psacalc delta numchamber, rmax(.10595)

qui reg use numchamber nummask i.($demographic base) ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base)
psacalc delta numchamber, rmax(.1092)

qui reg use numchamber nummask i.($demographic base) ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, vce(robust)
psacalc delta numchamber, rmax(.084)

eststo clear

********************************************************************************
*								WILD CLUSTER BOOTSTRAP					   *
********************************************************************************

eststo clear
	
eststo: qui reg use numchamber nummask, cluster(base)
boottest numchamber=0, boot(wild) weight(webb) nograph seed(8675309) level(90)
boottest nummask=0, boot(wild) weight(webb) nograph seed(8675309) level(90)

eststo: qui reg use numchamber nummask i.($demographic), cluster(base)
	estadd local Demog "Yes" 
boottest numchamber=0, boot(wild) weight(webb) nograph seed(8675309) level(90)
boottest nummask=0, boot(wild) weight(webb) nograph seed(8675309) level(90)

eststo: qui reg use numchamber nummask ally_relations i.($demographic), cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
boottest numchamber=0, boot(wild) weight(webb) nograph seed(8675309) level(90)
boottest nummask=0, boot(wild) weight(webb) nograph seed(8675309) level(90)

eststo: qui reg use numchamber nummask ally_relations orientation_officers i.($demographic), cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
boottest numchamber=0, boot(wild) weight(webb) nograph seed(8675309) level(90)
boottest nummask=0, boot(wild) weight(webb) nograph seed(8675309) level(90)

eststo: qui reg use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest i.($demographic), cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
boottest numchamber=0, boot(wild) weight(webb) nograph seed(8675309) level(90)
boottest nummask=0, boot(wild) weight(webb) nograph seed(8675309) level(90)

eststo: qui reg use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact i.($demographic), cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
boottest numchamber=0, boot(wild) weight(webb) nograph seed(8675309) level(90)
boottest nummask=0, boot(wild) weight(webb) nograph seed(8675309) level(90)

eststo: qui reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact i.($demographic), cluster(base) abs(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
boottest numchamber=0, boot(wild) weight(webb) nograph seed(8675309) level(90)
boottest nummask=0, boot(wild) weight(webb) nograph seed(8675309) level(90)

esttab using "${result}/CW_Wild.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic cov Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{15}{c}} \hline & \\ & & \multicolumn{13}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-15} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1 2 3 4 5 6 7) nomtit replace

eststo clear

********************************************************************************
*								MULTINOMIAL LOGIT							   *
********************************************************************************

eststo clear

eststo: mlogit use3 numchamber nummask, cluster(base) r base(0) rrr

eststo: mlogit use3 numchamber nummask 0.age 1.age 2.age 3.age 4.age 5.age 6.age 7.age 8.age 9.age 10.age 1.monthsoversea 2.monthsoversea 3.monthsoversea 4.monthsoversea 5.monthsoversea 6.monthsoversea 7.monthsoversea 8.monthsoversea 9.monthsoversea i.(school rankgrade), cluster(base) base(0) rrr
	estadd local Demog "Yes"

eststo: mlogit use3 numchamber nummask 0.age 1.age 2.age 3.age 4.age 5.age 6.age 7.age 8.age 9.age 10.age 1.monthsoversea 2.monthsoversea 3.monthsoversea 4.monthsoversea 5.monthsoversea 6.monthsoversea 7.monthsoversea 8.monthsoversea 9.monthsoversea i.(school rankgrade) ally_relations, cluster(base) base(0) rrr
	estadd local Demog "Yes"
	estadd local ForRel "Yes"

eststo: mlogit use3 numchamber nummask 0.age 1.age 2.age 3.age 4.age 5.age 6.age 7.age 8.age 9.age 10.age 1.monthsoversea 2.monthsoversea 3.monthsoversea 4.monthsoversea 5.monthsoversea 6.monthsoversea 7.monthsoversea 8.monthsoversea 9.monthsoversea i.(school rankgrade) ally_relations orientation_officers, cluster(base) base(0) rrr
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"

eststo: mlogit use3 numchamber nummask 0.age 1.age 2.age 3.age 4.age 5.age 6.age 7.age 8.age 9.age 10.age 1.monthsoversea 2.monthsoversea 3.monthsoversea 4.monthsoversea 5.monthsoversea 6.monthsoversea 7.monthsoversea 8.monthsoversea 9.monthsoversea i.(school rankgrade) ally_relations orientation_officers infocenter orient_meet war_interest, cluster(base) base(0) rrr
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"

eststo: mlogit use3 numchamber nummask 0.age 1.age 2.age 3.age 4.age 5.age 6.age 7.age 8.age 9.age 10.age 1.monthsoversea 2.monthsoversea 3.monthsoversea 4.monthsoversea 5.monthsoversea 6.monthsoversea 7.monthsoversea 8.monthsoversea 9.monthsoversea i.(school rankgrade) ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base) base(0) rrr
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"

eststo: mlogit use3 numchamber nummask 0.age 1.age 2.age 3.age 4.age 5.age 6.age 7.age 8.age 9.age 10.age 1.monthsoversea 2.monthsoversea 3.monthsoversea 4.monthsoversea 5.monthsoversea 6.monthsoversea 7.monthsoversea 8.monthsoversea 9.monthsoversea i.(school rankgrade base) ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base) base(0) rrr
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	
esttab using "${result}/CW_Multinomial.tex", style(tex) noomitted b(3) se(3) nonotes eform keep(numchamber nummask) stats(N aic cov Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{16}{c}} \hline & \\ & & \multicolumn{14}{c}{\textbf{Support for Second-Use of Chemical Weapons Against Japan (=1)}} \\ \cline{3-16} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1 2 3 4 5 6 7) nomtit replace

eststo clear

********************************************************************************
*									ADDED CONTROLS							   *
********************************************************************************

eststo clear

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	
eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact [pw=demweight], cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local DemWeight "Yes"
	
eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact training newmask, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local GasTrain "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact highprotect highconfid, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local Confid "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact correct, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local History "Yes"
	
eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact orient_talk, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local Cohesion "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact maskimpt chamberimpt gas_training, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local Important "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact physical_training military_training, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local MilTrain "Yes"
	
esttab using "${result}/CW_Controls.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic cov Demog ForRel Officer Info Local Base DemWeight GasTrain Confid History Cohesion Important MilTrain, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE" "\hspace{3mm}Demographic Weights" "\hspace{3mm}Chemical Training Details" "\hspace{3mm}Chemical Readiness" "\hspace{3mm}History of Use" "\hspace{3mm}Unit Cohesion" "\hspace{3mm}Chemical Training Importance" "\hspace{3mm}Military Training Importance")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{9}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-11} & \\ & & \textbf{Baseline} & & \multicolumn{7}{c}{\textbf{Additional Controls}}\\ \cline{2-3} \cline{5-11} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1 2) nomtit replace

eststo clear

********************************************************************************
*										PROBIT								   *
********************************************************************************

eststo clear
	
eststo: probit use numchamber nummask, cluster(base)

eststo: probit use numchamber nummask i.($demographic), cluster(base)
	estadd local Demog "Yes" 

eststo: probit use numchamber nummask ally_relations i.($demographic), cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"

eststo: probit use numchamber nummask ally_relations orientation_officers i.($demographic), cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"

eststo: probit use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest i.($demographic), cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"

eststo: probit use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact i.($demographic), cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"

eststo: probit use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact i.($demographic base), cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"

gen estimation_sample=(e(sample))
	
esttab using "${result}/CW_Probit.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic cov Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) nomtit replace

eststo clear
	
eststo: reghdfe use numchamber nummask if estimation_sample==1, cluster(base) noabs

eststo: reghdfe use numchamber nummask if estimation_sample==1, cluster(base) abs($demographic)
	estadd local Demog "Yes" 

eststo: reghdfe use numchamber nummask ally_relations if estimation_sample==1, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers if estimation_sample==1, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest if estimation_sample==1, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact if estimation_sample==1, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact if estimation_sample==1, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	
esttab using "${result}/CW_Main_Sample.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic cov Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) nomtit replace

eststo clear

********************************************************************************
*								COMPLIANCE WEIGHT							   *
********************************************************************************

eststo clear
	
eststo: reghdfe use numchamber nummask [pw=pr_complier], cluster(base) noabs
	estadd local Complier "Yes"

eststo: reghdfe use numchamber nummask [pw=pr_complier], cluster(base) abs($demographic)
	estadd local Demog "Yes" 
	estadd local Complier "Yes"

eststo: reghdfe use numchamber nummask ally_relations [pw=pr_complier], cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Complier "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers [pw=pr_complier], cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Complier "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest [pw=pr_complier], cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Complier "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact [pw=pr_complier], cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Complier "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact [pw=pr_complier], cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local Complier "Yes"
	
esttab using "${result}/CW_Comply.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic cov Complier Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}Pr(Compliance = 1)" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) nomtit replace

********************************************************************************
*										IPTW								   *
********************************************************************************

eststo clear
	
eststo: reghdfe use numchamber nummask [pw=iptw], cluster(base) noabs
	estadd local IPTW "Yes"

eststo: reghdfe use numchamber nummask [pw=iptw], cluster(base) abs($demographic)
	estadd local Demog "Yes" 
	estadd local IPTW "Yes"

eststo: reghdfe use numchamber nummask ally_relations [pw=iptw], cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local IPTW "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers [pw=iptw], cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local IPTW "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest [pw=iptw], cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local IPTW "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact [pw=iptw], cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local IPTW "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact [pw=iptw], cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local IPTW "Yes"
	
esttab using "${result}/CW_IPTW.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic cov IPTW Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}IPTW" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) nomtit replace

********************************************************************************
*										CEM									   *
********************************************************************************

eststo clear
	
eststo: reghdfe use numchamber nummask [pw=cem_weights1] if cem_matched1, cluster(base) noabs
	estadd local CEM "Yes"

eststo: reghdfe use numchamber nummask [pw=cem_weights1] if cem_matched1, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes" 

eststo: reghdfe use numchamber nummask ally_relations [pw=cem_weights1] if cem_matched1, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers [pw=cem_weights1] if cem_matched1, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest [pw=cem_weights1] if cem_matched1, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact [pw=cem_weights1] if cem_matched1, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact [pw=cem_weights1] if cem_matched1, cluster(base) abs($demographic base)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	
esttab using "${result}/CW_CEM1.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic, fmt(0 0 3 3) labels("Observations" "AIC")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ \textbf{Panel A:} & & \multicolumn{7}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) nomtit replace

eststo clear
	
eststo: reghdfe use numchamber nummask [pw=cem_weights2] if cem_matched2, cluster(base) noabs
	estadd local CEM "Yes"

eststo: reghdfe use numchamber nummask [pw=cem_weights2] if cem_matched2, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes" 

eststo: reghdfe use numchamber nummask ally_relations [pw=cem_weights2] if cem_matched2, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers [pw=cem_weights2] if cem_matched2, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest [pw=cem_weights2] if cem_matched2, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact [pw=cem_weights2] if cem_matched2, cluster(base) abs($demographic)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact [pw=cem_weights2] if cem_matched2, cluster(base) abs($demographic base)
	estadd local CEM "Yes"
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	
esttab using "${result}/CW_CEM2.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic cov CEM Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}CEM" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ \textbf{Panel B:} & & \multicolumn{7}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) nomtit replace

********************************************************************************
*										CMP									   *
********************************************************************************

eststo clear

cmp setup

eststo: cmp (use_comply= numchamber nummask) (complier = i.($demographic base) training orientation_officers), ind($cmp_cont $cmp_cont) cluster(base)

eststo: cmp (use_comply= numchamber nummask i.($demographic)) (complier = i.($demographic base) training orientation_officers), ind($cmp_cont $cmp_cont) cluster(base)
	estadd local Demog "Yes" 

eststo: cmp (use_comply= numchamber nummask i.($demographic) ally_relations) (complier = i.($demographic base) training orientation_officers), ind($cmp_cont $cmp_cont) cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"

eststo: cmp (use_comply= numchamber nummask i.($demographic) ally_relations orientation_officers) (complier = i.($demographic base) training orientation_officers), ind($cmp_cont $cmp_cont) cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"

eststo: cmp (use_comply= numchamber nummask i.($demographic) ally_relations orientation_officers infocenter orient_meet war_interest) (complier = i.($demographic base) training orientation_officers), ind($cmp_cont $cmp_cont) cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"

eststo: cmp (use_comply= numchamber nummask i.($demographic) nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact) (complier = i.($demographic base) training orientation_officers), ind($cmp_cont $cmp_cont) cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"

eststo: cmp (use_comply= numchamber nummask i.($demographic base) nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact) (complier = i.($demographic base) training orientation_officers), ind($cmp_cont $cmp_cont) cluster(base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	
esttab using "${result}/CW_CMP.tex", style(tex) b(3) se(3) nonotes keep(numchamber nummask) stats(N aic cov Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" nummask "& \\ Gas Mask Training") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{9}{c}} \hline & \\ & & \multicolumn{7}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-9} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1) nomtit replace

eststo clear

********************************************************************************
*									PEARL HARBOR							   *
********************************************************************************

gen pearlharbor=(monthsoversea>=8)

gen numchamberx=numchamber*pearlharbor
gen nummaskx=nummask*pearlharbor

eststo clear

eststo: reghdfe use numchamber nummask ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber numchamberx nummask nummaskx, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber numchamberx nummask nummaskx ally_relations, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber numchamberx nummask nummaskx ally_relations orientation_officers, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber numchamberx nummask nummaskx ally_relations orientation_officers infocenter orient_meet war_interest, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber numchamberx nummask nummaskx ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base) abs($demographic)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber numchamberx nummask nummaskx ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, cluster(base) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local Cluster "Yes"

eststo: reghdfe use numchamber numchamberx nummask nummaskx ally_relations orientation_officers infocenter orient_meet war_interest honolulu_contact, vce(robust) abs($demographic base)
	estadd local Demog "Yes"
	estadd local ForRel "Yes"
	estadd local Officer "Yes"
	estadd local Info "Yes"
	estadd local Local "Yes"
	estadd local Base "Yes"
	estadd local Cluster "No"
	
esttab using "${result}/CW_Combat.tex", style(tex) b(3) se(3) nonotes keep(numchamber numchamberx nummask nummaskx) stats(N aic cov Cluster Demog ForRel Officer Info Local Base, fmt(0 0 3 3) labels("Observations" "AIC" "\hline &\\ \textsc{Parameters}" "\hspace{3mm}Airbase Clustered SEs" "\hspace{3mm}Demographics" "\hspace{3mm}Postwar Foreign Policy" "\hspace{3mm}Officers' Leadership" "\hspace{3mm}Information Access" "\hspace{3mm}Local Contact" "\hspace{3mm}Airbase FE")) varlabels(numchamber "& \\ Gas Exposure" numchamberx "& \\ Gas Exposure x Pearl Harbor" nummask "& \\ Gas Mask Training" nummaskx "& \\ Gas Mask Training x Pearl Harbor") starlevels(* .10 ** .05 *** .01) nolines prehead(\begin{tabular}{l*{11}{c}} \hline & \\ & & \multicolumn{9}{c}{\textbf{Support for Using Chemical Weapons Against Japan (=1)}} \\ \cline{3-11} & \\ & & \textbf{Baseline} & & \multicolumn{7}{c}{\textbf{Additional Controls}}\\ \cline{2-3} \cline{5-11} & \\) posthead(\hline) prefoot() postfoot(& \\ \hline \end{tabular}) extracols(1 2) nomtit replace

eststo clear

********************************************************************************

clear
